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We present a derivation of a recently proposed theory for the time dependence of density fluctu¬ 
ations in stationary states of strongly interacting, athermal, self-propelled particles. The derivation 
consists of two steps. First, we start from the equation of motion for the joint distribution of parti¬ 
cles’ positions and self-propulsions and we integrate out the self-propulsions. In this way we derive 
an approximate, many-particle evolution equation for the probability distribution of the particles’ 
positions. Second, we use this evolution equation to describe the time dependence of steady-state 
density correlations. We derive a memory function representation of the density correlation function 
and then we use a factorization approximation to obtain an approximate expression for the memory 
function. In the final equation of motion for the density correlation function the non-equilibrium 
character of the active system manifests itself through the presence of a new steady-state correla¬ 
tion function that quantifies spatial correlations of the velocities of the particles. This correlation 
function enters into the frequency term, and thus it describes the dependence of the short-time 
dynamics on the properties of the self-propulsions. More importantly, the correlation function of 
particles’ velocities enters into the vertex of the memory function and through the vertex it modifies 
the long-time glassy dynamics. 

PACS numbers: 82.70.Dd, 64.70.pv, 64.70.Q-, 47.57.-s 


I. INTRODUCTION 

There is a lot of interest in the static and dynamic 
properties of active matter systems EH1]- These sys¬ 
tems consist of particles that are self-propelled, and thus 
can move autonomously. The reason for this interest is 
twofold. First, there exists a great variety of biological 
and physical systems with self-propelled constituents. In 
this work we consider only a subset of these systems con¬ 
sisting of assemblies of particles without aligning interac¬ 
tions. Biological examples of such systems include some 
motile bacteria and physical systems are represented by 
synthetic Janus colloidal particles. The second reason 
for the interest in active matter systems is that they pro¬ 
vide fascinating examples of intrinsically non-equilibrium 
driven systems with very rich and unusual phenomenol¬ 
ogy. From the point of view of a theoretical analysis 
active systems are a little easier to study than sheared 
systems, which used to be a workhorse of non-equilibrium 
statistical mechanics Q- This follows from the fact that 
in active matter systems the driving is applied locally 
and, on average, isotropically. In contrast, sheared sys¬ 
tems are intrinsically anisotropic. 

Recently, it has been realized that dense active systems 
can exhibit glassy dynamics. The analogy between the 
dynamics in dense cell layers and the dynamics in a fluid 
approaching a glass transition was noted by Angelini et 
al. [5}. This study inspired a simulational investigation 
of a two-dimensional dense active system with aligning 
interactions in which a jammed phase was identified 0. 
Glassy dynamics and the active glass transition was ana¬ 
lyzed in a more general setting by Berthier and Kurchan 
J]. They studied dynamic behavior of a simple model ac¬ 
tive system inspired by the spherical p-spin model. Their 


investigation showed that activity can shift the onset of 
glassy behavior but it does not always destroy the glass 
transition, unlike for similar models of sheared systems 
[§]. Berthier and Kurchan’s theoretical analysis was fol¬ 
lowed by several computational studies of strongly inter¬ 
acting systems of self-propelled particles without aligning 
interactions, which demonstrated glassy behavior of ac¬ 
tive matter. First, Ni et al. 0 simulated a system of 
active Brownian hard spheres and found that, with in¬ 
creasing magnitude of the self-propulsion, the long-time 
dynamics speeds up whereas the local structure, moni¬ 
tored through the steady-state structure factor, becomes 
less pronounced. Next, Berthier EH simulated a slightly 
different, fully athermal (i.e. devoid of thermal Brownian 
motion) system of active hard disks and found that, with 
increasing departure from equilibrium, at intermediate 
volume fractions the dynamics could be non-monotonic 
but, again, the apparent glass transition is shifted to¬ 
wards higher volume fractions. Wysocki et al. EH sim¬ 
ulated a dense system of active Brownian particles with 
continuous interactions. More recently, Mandal et al. 
EH simulated the classic model glassy system, the Kob- 
Andersen Lennard-Jones binary mixture E3> i 11 which 
one of the components was endowed with self-propulsion. 
In qualitative agreement with Ni et al., they found that, 
upon increasing the magnitude of the self-propulsion, the 
long-time dynamics speeds up and that, beyond a crit¬ 
ical strength of the self-propulsion, the glass phase dis¬ 
appears. Finally, Fily et al. EH investigated the phase 
diagram of active Brownian harmonic spheres as a func¬ 
tion of density, activity and noise, and identified a glassy 
phase in the high density, small self-propulsion speed 
regime. 

The simulational investigations of active glassy dy- 
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namics stimulated interest in theoretical descriptions of 
the dynamics of dense systems of self-propelled particles. 
The main goal of these descriptions is to understand the 
effect of the non-equilibrium character of the active sys¬ 
tem on the glassy dynamics and the glass transition. We 
recall that, in spite of an extensive theoretical effort, 
a microscopic and quantitatively correct description of 
glassy dynamics of thermal “passive” systems remains a 
challenge [Tsj . Thus, we should not expect to obtain a 
microscopic and quantitatively accurate theory for much 
more complicated active systems. However, we might be 
able to get some insight into the role of activity. The 
precedent for this is the application of the well-known, 
albeit far from perfect, theory of the glassy dynamics, 
the mode-coupling theory [lf|, to the description of col¬ 
loidal systems with short-range attractive interactions 
E2- The theory, while in general not quantitatively cor¬ 
rect, was able to account for the effect of attractions and 
to predict a re-entrant glass transition. Its qualitative 
predictions were subsequently confirmed by experiments 
jl8| . We note that proposed theoretical descriptions of 
the dynamics of dense active systems use approaches sim¬ 
ilar the mode-coupling theory. This is hardly surprising 
since the mode-coupling theory provides one of very few 
microscopic descriptions of the glassy dynamics in ther¬ 
mal systems [n|. In addition, the virtues and flaws of 
this theory are quite well understood, which allows one to 
focus on phenomena that the theory describes reasonably 
well. 

The first microscopic theory for the glassy dynamics 
of dense systems consisting of spherically symmetric self- 
propelled particles with non-ali gnin g interactions was put 
forward by Farage and Brader [201 ]. This theory is a gen¬ 
eralization of an earlier mode-coupling theory for sheared 
glassy colloidal fluids proposed by Fuchs and Cates [U- 
Farage and Brader considered a system of self-propelled 
particles that were also subjected to thermal noise. They 
used the integration-through-transients approach intro¬ 
duced by Fuchs and Cates. In this approach one assumes 
that in the infinitely distant past the system was in an 
equilibrium state and then the drive, in this case the ac¬ 
tivity, was turned on. Using this procedure one can, in 
principle, calculate both equal time and dynamic proper¬ 
ties of active systems. However, it is not clear how to use 
this approach to describe a fully athermal system which 
does not have an equilibrium state without the drive. 

Recently, we proposed a different microscopic theory 
for dense athermal active systems [|I| ■ Like the theory of 
Farage and Brader, our theory relies upon a factorization 
approximation, and thus it falls into a general category of 
mode-coupling-like approaches. There are, however, sig¬ 
nificant differences between our theory and that of Farage 
and Brader. First, our approach incorporates the influ¬ 
ence of correlations between velocities of different parti¬ 
cles on the dynamics of the active system. This effect is 
neglected in the theory of Farage and Brader. Second, 
we use a projection operator defined through a steady- 
state distribution function. This approach avoids using 


the integration-through-transients procedure and allows 
us to describe dynamics of fully athermal systems. 

In Ref. (2lJ we presented an outline of the derivation 
of our theory. Here we present the details of the deriva¬ 
tion and discuss all the approximations used. To make 
the present contribution self-contained we include some 
of the discussion from the appendices of Ref. [21]. In 
addition, we present a derivation of the theory for the 
time dependence of the tracer density fluctuations. 

The main object of our theory is the (collective) inter¬ 
mediate scattering function, F(q;t), which describes the 
time dependence of the collective density fluctuations, 

1 ^ e - iq ' riW E eiq rj(0) ^ ' C 1 ) 

Here and in the following the brackets (...) denote av¬ 
eraging over a steady-state distribution of positions and 
self-propulsions. We will also consider a closely related 
function, the self-intermediate scattering function, which 
describes the time dependence of the tagged particle 
(tracer) density fluctuations, F s (q;t), 

F s (q;t) = 1 ^ e - iq -(-(*)- r ‘(°))^ . (2) 

The latter function allows us to calculate quantities that 
are usually monitored in computer simulations,! e the 
relaxation time t q defined through the relation [22| 

F a {q\T a ) = e~ x , (3) 

and the (long-time) self-diffusion coefficient D defined 
through the hydrodynamic limit of the self-intermediate 
scattering function, 

lim F s (q;t) ~ exp (- Dq 2 t ) . (4) 

t—^oo,q^O,q 2 t—const. 

The paper is organized as follows. In Section [TT] we 
present and motivate the model active system that we 
analyze in the reminder of the paper. In Section Mil we 
discuss the main assumption on which our approach re¬ 
lies, the absence of steady-state currents. In Section IIVI 
we present the derivation of the approximate equation 
of motion for the many-body distribution of the posi¬ 
tions of the active particles. This is followed by a brief 
discussion (Section 0 of the short-time behavior of the 
scattering functions which illustrates the importance of 
the correlations between self-propulsions and positions. 
The derivations of the memory function representation 
and of the approximate expression for the memory func¬ 
tion are presented in sections fVTl and lYlIl respectively. In 
Section IVIIII we discuss general conclusions that can be 
drawn from the final equation for the time dependence 
of the density correlations and in Section IIXI we summa¬ 
rize our findings and outline directions of future work. In 
the first appendix, Appendix E we discuss an alterna¬ 
tive equation of motion for the many-body distribution 





3 


of the particles’ positions. In the second appendix, Ap¬ 
pendix Q3] we present the derivation of the theory for the 
time dependence of the tagged particle (tracer) density 
fluctuations. 


II. MODEL ACTIVE SYSTEM 

We consider a system of N interacting, self-propelled 
particles in a volume V. The average density is p = N/V. 
The particles interact via a spherically symmetric poten¬ 
tial V (r). They move in a viscous medium that is char¬ 
acterized by the friction coefficient of a single particle, 
which we denote by £o- We assume that the friction felt 
by a particle is independent of the particle density and 
configuration, and thus we neglect hydrodynamic inter¬ 
actions [23} • Each particle moves under the combined 
influence of the interparticle force derived from the po¬ 
tential V (r) and a self-propulsion force [24j • 

In the active Brownian particles model in Refs. j25l.(26j 
it is assumed that the magnitude of the self-propulsion 
force is constant whereas its direction changes via rota¬ 
tional Brownian motion. If needed, to distinguish this 
model from other models of active particles we will refer 
to it as the rotational diffusion active Brownian particles 
model. 

We use a slightly different model of active particles, in¬ 
troduced in Ref. [27| and then analyzed in Ref. [21] . An 
essentially identical model was independently introduced 
and studied by Maggi et al. (28 ]. 

Specifically, we assume that the self-propulsion force 
does not have a constant magnitude and it evolves in 
time according to the Ornstein-Uhlenbeck process. We 
note that the rotational diffusion active Brownian parti¬ 
cles model and our model can be considered to be limit¬ 
ing cases of a more general model in which the direc¬ 
tion of the self-propulsion force randomly rotates and 
its magnitude fluctuates. Most likely, neither limiting 
case perfectly describes active particles studied in exper¬ 
iments. The advantage of our model is that the equa¬ 
tions of motion for the self-propulsion are linear, which 
makes some theoretical considerations somewhat easier. 
In particular, we showed that the complete steady-state 
distribution can be found for a single particle in a har¬ 
monic potential for our model [23] - We note that our 
model can be considered to be a continuous time version 
of the model used by Berthier 0- Finally, we mention 
that a mapping was derived between the rotational dif¬ 
fusion active Brownian particles model and our model 
[2dl [29[ | . This mapping preserves the self-propulsion au¬ 
tocorrelation function. However, higher moments of the 
self-propulsion are different. 

We assume that the self-propelled particles are large 
enough so that any Brownian motion due to the ther¬ 
mal fluctuations of the viscous medium (solvent) can be 
neglected. Thus, the system is purely athermal and the 
particles move under the sole influence of the interparti¬ 
cle interactions and self-propulsion. 


The above qualitative discussion translates into the 
following equations of motion for the positions and self¬ 
propulsions, 

D = £ 0 -1 [fi+Fi], (5) 

fi = -T p _1 fi +Vi- (6) 

In Eq. ©, r, is the position of particle i, £0 the friction 
coefficient of an isolated particle, f \ is the self-propulsion 
acting on particle i and F i is the force acting on particle 
i originating from the interactions, 

= ( 7 ) 

In Eq. ([6]) , t p is the persistence time of the self-propulsion 
and r] i is an internal Gaussian noise with zero mean and 
variance 

(‘ni(t) T lj(t ')) noise = 2 D f I6 i:j 6(t - t'), (8) 

where (•■■} noiso denotes averaging over the noise distri¬ 
bution, Df is the noise strength and I is the unit ten¬ 
sor. Without interactions, particles evolving according 
to Eqs. m perform a persistent random walk with the 
mean-square displacement |30) 

((r i(t) ~ D(0)) 2 ^} = (t + T p (e“ t/rp - 1)) . (9) 

According to Eq. ([9]), the short-time motion is ballistic, 

— r i(0)) 2 \ ~ 3-^2~^f 2 (10) 

' 'so 

and the long-time motion, l>r p , is diffusive with diffu¬ 
sion coefficient Do, 

Do = D f Tlie o . ( 11 ) 

Comparing expression m with the well-known for¬ 
mula for the diffusion coefficient of a Brownian particle 
moving in a viscous medium with friction constant £ 0 , 
^Brownian = T /£o (we use units such that the Boltzmann 
constant ks = 1), we can define the single-particle effec¬ 
tive temperature [23 ] . 

T e S = A)£o = DfTp/t; o- (12) 

The single particle effective temperature could be used 
as one of the independent control parameters, together 
with the number density p and the persistence time r p . 

Eqs. ([5][G]) are a convenient starting point for com¬ 
puter simulations. In theoretical considerations it is more 
useful to describe the system using probability distri¬ 
butions and their associated evolution equations. The 
most fundamental description of our system, equivalent 
to equations of motion m , uses the joint iV-particle 
probability distribution of positions and self-propulsions, 
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Pjv(ri, fi,rjv, fjv! t). This distribution evolves in time In Eq. (1201) n(q) is the Fourier transform of the micro- 
with evolution operator fl scopic density, 


dtPjv-(ri,fi,...,rjv,fjv;£) = ClP N (ri, fi, r N , fjv; t), 

_ (13) 

which can be derived [30| from equations of motion ([5][6]), 

SI = £ V, • (f, + F.) + £ ±. (r-'f, + D f ±). 

i i ' ' 

(14) 

We recall that for non-interacting active particles Eqs. 
mm are formally equivalent to the Fokker-Planck 
equation that describes the motion of non-interacting 
Brownian particles on a time scale on which their ve¬ 
locity relaxation can be observed [3Cj. Thus, a system of 
non-interacting active particles is formally equivalent to 
particles moving under the influence of thermal solvent 
fluctuations. This equivalence is absent for a system of 
interacting active particles. In particular, it should be 
emphasized that according to Eqs. m or, equivalently, 
Eqs. (I13H14I) . self-propulsions evolves autonomously, in¬ 
dependently of the configuration of the particles. 

We will assume that our system can achieve a station¬ 
ary state. In other words, there exists a steady-state 
probability distribution Pjy (iq, fi,..., rjv, fjv) such that 

^P/v( r ii fi, rjv, fjv) = 0. (15) 

We emphasize that, in general, the joint steady-state 
distribution of positions and self-propulsions does not 
factorize into a product of steady-state distributions of 
particle positions and self-propulsions, 

P s N s (ri, fa,rjv, f N ) * Pn( n,.... r Jv)Pjv (fi>fjv), 

(16) 

where P^(ri,..., rjv) and Pjy (fi,..., fjv) are the steady- 
state distributions of positions and self-propulsions, 

Pjv ( r i, Uv) = J fi,..., r N , fjv),(17) 

P.v(f|,-,fjv) = J dr 1 ...dr N P^(r 1 ,ii, ..., rjv, fjv)(18) 

In general, neither the joint steady-state distribution 
Pjv(ri, fi,..., r^r, fjv) nor the steady-state distributions of 
positions Pf$ (ri, ..., rjv) are known exactly (for approx¬ 
imate theories for the latter distribution see Refs. (Hj, 
l29jl. However, the steady-state distribution of self¬ 
propulsions has a simple form, 

JTO.-.ftv) oce xp(-§^)- (19) 

Evolution operator d allows us to rewrite the defi¬ 
nitions of the intermediate scattering functions m , 

F{q\ t) = ^ (n(q) exp (fit) n(-q)), (20) 

(21) 


n(q)=5>-^ r ‘, (22) 

i 

and in Eq. (1211) n s (q) is the Fourier transform of the 
microscopic tagged particle (tracer) density, 

n s (q)=e- iqri . (23) 

Note that we chose the tagged particle to be the particle 
number one; this choice is made arbitrarily and for conve¬ 
nience only. In Eqs. (12011211) H is the evolution operator 
m, and brackets (...) denote averaging over the joint 
steady-state distribution of positions and self-propulsions 
m- We emphasize that in Eqs. (12011211) and in all sim¬ 
ilar formulas the steady-state distribution stands to the 
right of the quantity being averaged, and all operators 
act on it too. 


III. THE MAIN ASSUMPTION: ABSENCE OF 
CURRENTS 


The fundamental difficulty posed by driven systems 
originates from the absence of detailed balance. The lack 
of detailed balance allows for the existence of non-trivial 
currents. It has been argued that at a mesoscopic level, 
in some systems without aligning interactions, i.e. after 
some coarse-graining (i.e. above a certain length and 
time scale), there are no currents and the active system 
becomes equivalent to a passive system [31.]. 

Our main assumption is that in our model system, in 
the steady state the currents vanish after integrating out 
the self-propulsions. This assumption, together with ad¬ 
ditional coarse-graining over time, will allow us to ap¬ 
proximate our system by a passive system with detailed 
balance. 

A similar assumption was implicitly used by Farage 
and Brader. However, as discussed in the next section, 
their theory and ours use different passive systems to 
approximate an active system. In particular, we retain 
correlations between (overdamped) velocities of different 
particles. 

To make our assumption explicit we first rewrite the 
equation of motion for the joint probability distribution 
of positions and self-propulsions, Eq. m, in the form 
of a continuity equation, 


diPjv(ri,fi,...,rjv,fjv;t) = 

Vj • jj(ri,fi,...,r^r,fjv; t) 

i 

•ji(ri,fi,-,rjv,fjy;£), 


l 

where current densities are defined as 


(24) 


ji(ri,fj, ■■■,rjv,fjv; t) = 

CcT 1 (Fi + fi) Pv(ri, fi,..., rjv, fjv; t), 


F s (q; t) = (n s (q) exp (fit) n s (- q)). 


(25) 
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jf ( r i,fi, ...,r N ,f N ;t) = (26) 

- + D f-^j Pjv( r i> fi, •••, r N ,f N -, t). 

Current densities (I25H26I) are microscopic quantities, 
which may be non-zero in a system without detailed bal¬ 
ance. Our main assumption is that in the steady state, 
the current density in the position space, integrated over 
self-propulsions, vanishes, 

jf(ri,-,r N ;t) = Co J dfi-dfjvjf (ri,fi, ...,r N ,f N ;t) 

= Co -1 J dh.-.dfN [F, + fi] P^ s (n, f 1; r N , f N ) 

= 0. (27) 

Assumption m implies that the local steady-state av¬ 
erage of the self-propulsion is equal to the negative of the 
force, 

<fi>i ss = -F<, (28) 

where the local steady-state average is defined as 

<••■>!»= (29) 

Eq. (17S1) could be interpreted as stating a balance of the 
self-propulsion of particle i and the total potential force 
acting on this particle, for a given configuration of the 
system. 

We are not aware of any study that specifically fo¬ 
cused on the existence of non-trivial steady-state currents 
in high density active systems without aligning interac¬ 
tions jH]. We note that the assumption (1271) is made at 
the level of V-particle quantities. Thus, its direct sim- 
ulational verification seems rather difficult. However, it 
might be possible to define and measure reduced (few- 
particle) current densities. Alternatively, it might be 
possible to analyze theoretically and then verify com¬ 
putationally the consequences of the presence/absence of 
currents. For example, one could try to derive a theory 
for an effective temperature of the active system and then 
to compare it with a recent simulational investigation 
[32| . Finally, one could try to develop an approximate 
expression for the A-particle joint steady-state distribu¬ 
tion of positions and self-propulsions and investigate the 
validity of Eq. (1271) directly. Preliminary results from 
such a project Q suggest that, in general, steady-state 
currents exist and thus Eq. (EH) constitutes an approx¬ 
imation. We will return to the consequence of the as¬ 
sumption made in Eq. (1271 ) in future work. 


IV. EFFECTIVE EQUATION OF MOTION FOR 
THE DISTRIBUTION OF PARTICLE POSITIONS 


As emphasized in Sec. El in our model active sys¬ 
tem the evolution of the self-propulsions does not depend 


on the configuration and density of the particles. On 
the other hand, for any interparticle interactions with a 
strongly repulsive short-range part of the interaction po¬ 
tential, in a dense system the positions of the particles 
change slowly. Intuitively, upon increasing the strength 
of the interactions in a dense active system, one expects 
qualitatively similar slowing down as in a passive sys¬ 
tem. By and large, this expectation has been borne out 
by computer simulations [3, [jX], |l2|, 1EI albeit with some 
caveats. 

Thus, it seems reasonable to assume that in strongly 
interacting systems of active particles without aligning 
interactions the self-propulsions relax faster than the po¬ 
sitions of the particles. This separation of times scales 
suggests that it should be possible to derive an approxi¬ 
mate equation of motion for the probability of particles’ 
positions only. We expect that such an equation should 
become progressively more accurate in the strong inter¬ 
action limit. 

We should recall that our goal is to describe the dy¬ 
namics of density fluctuations in the steady state of our 
model active system. Thus, we only need an approxi¬ 
mate equation of motion for the probability of particles’ 
positions in the vicinity of the steady state. For this rea¬ 
son, in our considerations we introduce the steady-state 
distribution and, effectively, we use the gradient of the 
logarithm of this distribution with respect to the posi¬ 
tion of a given particle as a (normalized) effective force 
acting on this particle (see Eq. (©). This construction 
is different from the approach of Farage and Brader [20j] 
who concentrate on transient density fluctuations. 

To facilitate the derivation of the equation of motion 
for the probability of particles’ positions we introduce 
a projection operator that acts on an V-particle prob¬ 
ability distribution of self-propulsions and positions and 
projects it on a local steady-state distribution, i.e. on a 
distribution in which self-propulsions have a steady-state 
distribution for a given sample of positions, 


Vi ss P N (ri,fi, rjv, fjv; t) 

Pn (ri,fi,-,rjy,fjv) 

Pn( ri,-,rjv) 
P$(r 1 ,f l ,...,r N ,i N ) 
P|f(ri,...,rjv) 


j df 1 ...df N P N (r 1 ,f 1 ,...,r N ,f N \t) 
P JV (r 1 ,...,r JV ;t). (30) 


We note that by integrating Pi ss Piv(ri, fi,rjv, fjv; t) 
over self-propulsions we get the probability distribution 
of particles’ positions, P/v(ri,..., r^r; t). 

Next, we define the orthogonal projection, 


Qlss = I ~ Piss, (31) 


and write down equations of mo¬ 
tion for PissPjv(ri,fi,...,rjv,fjv;t) and 
QusPjv(ri,fi,..., rjv, fjv; t), 

5tPi ss Pjv(ri,fi,...,rjv,fjv;t) = 

Pi ss f2Pi ss Pjv(ri, fi,..., rjv, fjv; t) 
+Pi ss C2QissPiv(ri,fi,..., rjv, f N -1), (32) 
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<9iQi S s-P/v( r i; fij ■■■) r Ar, f/v; t) = 

Qi ss nV lss P N (ri,fi, rjv,fjv; t) 
+fii»fiQi»iV(ri,fi ) rjv, fjv;i)- (33) 

Since our goal is to calculate the intermediate scat¬ 
tering functions, Eqs. HD , which are functions of 

I 


positions only, we can assume that 

Qiss-Pzv(ri, fi,r/v, i N ; t = 0) = 0. (34) 

Then we can solve Eqs. (132113311 for the Laplace trans¬ 
form, CT, of d t V\s S Pisr(ri, fq,..., r N , fjvi t) and we obtain 


CT [9t'Pi ss Pjv(ir, fi,..., rjy, f^v; t)] (z) = 


V\ ss nP\ ss + Pi ss nQ ls 


i 


' z - QissOQi, 

The first term inside the brackets on right-hand-side of Eq. ( 1551 ) reads 

Pn (ri,fi,...,rjv,fjv) -_i 


-Qlss^Pls, 


Piss-Pjv(ri,fi, ...,rjv,fjv;z)(35) 


V ls JlVi ss P N (r-j, f]..... r jY , fv! z) = 


xV 


Pn (ri,...,rjy) 

Pjy(ri,...,rjy;z) 

P^(ri,...,riv) 


1 / dfi...dfiv [Fi + P£(n, f 1; ..., rjv , fjv) 


(36) 


We see that if current densities vanish in the steady state, Eq. m, this term vanishes. Furthermore, one can show 
that 


Qi ss nr lss P N (z) = E (fi - (fi) lss ) i^(ri, fi,..., rjv, fjv) 


P N (r 1 ,...,r N ;z) 


_ 1 P^(n,...,r N ) _ 


and 


P\ ss ClQ lss ... = E ' / df i-^ ( f * - < f ‘U - 


(37) 


(38) 


r 


So far, we have not made any approximations. To pro¬ 
ceed, we will need to deal with projected evolution op¬ 
erator Qi ss flQi ss in Eq. ( 1551 ) . This operator describes 
evolution in the space orthogonal to the local steady- 
state space. The simplest possible approximation is to 
assume that this evolution is entirely due to the free re¬ 
laxation of the self-propulsions. In this case Qi ss fIQi S s is 
approximated as follows 

Qis S ^Qi ss « E E L-if. + Df A) . (39) 

i— 1 1 \ l ' 

We note that approximation (1391) is equivalent to assum¬ 


ing that the time-scale of the self-propulsion’s relaxation 
is the same in non-interacting and interacting systems. 
This can be approximately valid for our model active 
system but, in particular, it would be an unreasonable 
approximation in the presence of aligning interactions. 
We recall that in our system the evolution of the self¬ 
propulsions is independent of the positions of the parti¬ 
cles. Thus, the approximation ( 1551 ) neglects the influence 
of the correlations between self-propulsions and positions 
on the evolution of the self-propulsions. Combining this 
approximation with Eqs. (13711381) we get the following 
approximate equality 


Vlss^Qlss (z - Slss^Qlss)- 1 Qi ss ^Pis S Pn(Vi, f 1, ..., V N , f N] Z ) « ^^ E V > ’ f dfl-.di N (f,; - ($>jJ 


-P|f( r i, r/v) 


N 


-En?T ( T p lf :+ D f 


dij 
j=i J 


Now, we expand 


<9f, 


E(ft-(ft)iJ^ s (ri,fi ! ...,r JV ,fiv) 


V; 


-P/v( r i, rjy; z) 
Pn( ri,...,rjv) 


( 40 ) 


* - Eili m- ( T i lf i + D fm) 


i -i 


and integrate by parts. Finally, we integrate both sides of the 
resulting equation over self-propulsions and get the following expression for the Laplace transform of 9 t Pjv(iq,..., r at; f) 

CT[d t P N (n,...,r N -t)\ (z) = C 0 - 2 E Vi • (z + Tp 1 ) 1 - (fi)^ (f # ) J ■ [-Ff + V,] P N (r 1 ,...,r N] z). (41) 


h3 
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where F° ff is the (normalized) effective force acting on 
particle j in the steady state, 

Ff = V J \nP s N s (r 1 ,...,r N ). (42) 

The right-hand-side of Eq. m defines the effective evo¬ 
lution operator n off (,G), 

^ eff W = 4' 2 E v *-(^ + r p -1 )" 1 

«f^)iss - &>!*, <f,'> J • [-Ff + V,] . (43) 

We recall that in approximation (l39l) the influence of the 
correlations between self-propulsions and positions on the 
evolution of the self-propulsions was neglected. In the ef¬ 
fective evolution operator fl51) the influence of the same 
correlations on the (much slower) evolution of the parti¬ 
cles’ positions is included through the steady-state cor¬ 
relations of the self-propulsions, (f)fj) lss - (fi) lss (fj) lss - 

Furthermore, we note that in the z — > 0 limit (which 
corresponds to coarse-graining in time) the evolution op¬ 
erator, Eq. (l43l) . becomes formally similar to the evo¬ 
lution operator for a system of Brownian particles with 
hydrodynamic interactions [23| . This is the formal ex¬ 
pression of the equivalence of a coarse-grained active sys¬ 
tem ( i.e. in the space of positions only and on long time 
scales) and a passive system. 

Farage and Brader [201 ] implicitly assumed the equiv¬ 
alence of a similarly coarse-grained active system {i.e. 
in the space of positions only and on long time scales) 
and a passive system without hydrodynamic interactions. 
Specifically, in order to obtain the effective evolution op¬ 
erator used in Ref. [2C| from our operator (H3l) . one needs 
to replace local steady-state correlations between the self¬ 
propulsions by the mean-square of the self-propulsion, 

( f * f j}i ss - ( f i)iss ( f j)i ss -»■ (f^noise % = D f T p S i:j I, where 
I is the unit tensor, and one needs to take the z —>■ 0 
limit. Thus, our theory and the theory of Farage and 
Brader use different passive systems to approximate a 
coarse grained active system. 

The most important assumption made in this section 
was approximating Qi ss ^Qiss by the free relaxation of the 
self-propulsions. We shall see later that this assumption 
allows us to recover the correct short time dependence of 
the intermediate scattering function. 

We recall once again that according to our model’s evo¬ 
lution equations, Eqs. (151161) . the self-propulsions evolve 
independently of the positions. However, in general, 
due to the coupled steady-state distribution of the self¬ 
propulsions and positions, the characteristic time for the 
transition from the short-time ballistic to diffusive mo¬ 
tion (which is equal to r p for non-interacting particles, 
see Eq. ([9])) depends on the strength of the interparticle 
interactions. We referred to this change in Ref. [2l| as a 
renormalization of the relaxation rate of the distribution 
of the self-propulsions by the interactions. As we will 
see later, approximating Qi ss HQi ss by the free relaxation 
of the self-propulsions is equivalent to neglecting of any 
such renormalization. 


To improve upon this approximation one could try the 
alternative approach presented in Appendix [A] This ap¬ 
proach will likely lead to a more accurate description of 
the short-time dynamics. However, this advantage is out- 
weighted by its computational complexity. In the remain¬ 
der of this paper we will use the approximate effective 
evolution operator (PT51) . 


V. SHORT-TIME DYNAMICS: IMPORTANCE 
OF CORRELATIONS BETWEEN 
SELF-PROPULSIONS AND POSITIONS 


Here we briefly discuss the short-time dynamics of the 
scattering functions. First, we derive and discuss the ex¬ 
act expressions for the second time derivatives of the scat¬ 
tering functions. We emphasize the importance of the 
correlations between the self-propulsions acting on the 
particles and their positions for the correct description 
of the short-time dynamics. Second, we point out that 
the effective evolution operator (T451) leads to the correct 
results for the second time derivatives of the scattering 
functions. 

To evaluate the short-time behavior of the intermediate 
scattering function we expand expression © in powers 
of t , 

(n(q)n(-q)) + (n(q)fin(-q)) 

+ ^N + - ( 44 ) 

The first term at the right-hand-side of Eq. © is the 
steady state static structure factor, 

-%) = Mq)n(-q)), (45) 

the second term vanishes due to symmetry, and the third 
term gives 

^(n(q)H 2 n(-q)) = 

= -~^~ UJ \\ (?)• ( 46 ) 

In Eq. (1461) . W||(g) is a function that quantifies corre¬ 
lations of velocities of individual particles, 


W||(g) = 
1 „ 




o 



(47) 

q 


The interpretation of function W||(g) comes from the fact 
that (fi + Fj) is the velocity of particle i for our over¬ 
damped system, see Eq. ©■ Computer simulation re¬ 
sults [2lJ show that, at constant single particle effective 




temperature T e g, Eq. m, with increasing the persis¬ 
tence time ( i.e . with increasing departure from equi¬ 
librium), W||(g) gradually develops significant oscillations 
while its overall magnitude decreases. 

To summarize, the short-time expansion reads, 

+ (48) 

This should be compared to the standard result for a 
colloidal suspension, 

F BD (q;t) = S(q)-q 2 tD 0 H(q) + ... (49) 

where Dq is the diffusion coefficient of a single isolated 
Brownian particle and H ( q ) is the hydrodynamic factor 
[2dt and the Newtonian dynamics result, 

F nd (,;() = S( 5 ) - ^21 + ... (50) 

2 TO 

where fee is the Boltzmann constant, T is the tempera¬ 
ture and to is the particle mass. 

The analogous calculation for the self-intermediate 
scattering function gives 

F s (q;t) = l-^y-W||(oo) + ... (51) 

whereas for a colloidal system one gets 

F BD (q;t) = l-q 2 tD 0 H(oo) + ... (52) 

where D 0 is the diffusion coefficient of a single isolated 
Brownian particle and F[{p o) = linig-^oo H(q), and the 
result for a Newtonian system reads 

Frw.t) = + ... ( 53 ) 

The above results show that our active system has fea¬ 
tures of both colloidal and Newtonian systems. Its short- 
time dynamics involves a non-trivial steady-state corre¬ 
lation function W||(g) quantifying spatial correlations of 
particles’ velocities which plays a role similar to that of 
the hydrodynamic factor H(q). On the other hand, the 
short-time dynamics of the active system is ballistic, like 
in a Newtonian system. However, the origin of this bal¬ 
listic behavior is self-propulsion rather than inertia. 

It can be showed that the same results for the short- 
time dynamics are obtained if one starts from expressions 
for the scattering functions in terms of the effective evo¬ 
lution operator fl51) . The difference between the exact 
expressions and those involving the effective evolution 
operator appears at t 3 order. 


VI. MEMORY FUNCTION REPRESENTATION 

In this section we rewrite the formal expression (l20l) 
for the intermediate scattering function in terms of the 
so-called frequency matrix and irreducible memory ma¬ 
trix. The latter quantity contains all the (unknown) non¬ 
trivial dynamic information about the system. The re¬ 
sulting expression for the density correlation function in 
terms of the frequency matrix and the memory matrix is 
known as the memory function representation. 

We start by re-writing of the Laplace transform of the 
intermediate scattering function in terms of the approx¬ 
imate evolution operator (T451) . 

CT[F(q-t)\ (z) = F(q; z) = 

(n{ q) (z - ny 1 n(-q)^ = 

(n{q) (z-n eS (z)) 1 n(-q)^ . (54) 

Here (...) r denotes averaging over the steady-state distri¬ 
bution of particles’ positions. 

To derive the memory function representation of 
F(q ; z) we use the projection operator approach HE HE 
|39(. We define a projection operator on the microscopic 
density 

V n = ... n(- q)) r (n(q)n(-q))" 1 (n( q)... . (55) 

We emphasize that projection operator V n is defined in 
terms of the steady-state distribution, unlike in the ap¬ 
proach of Farage and Brader j20|. Next, we use the iden¬ 
tity 


z-n eff (z) z-n eS (z)Q n v 

+ * - (z) Q n neS[z)Vn z - (z) ’ 

where Q n is the projection on the space orthogonal to 
that spanned by the microscopic density, 

Q n =l~V n , (57) 

to rewrite the Laplace transform of the time derivative 
of NF(q\t) in the following way 
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CT[d t NF(q-t)\(z) = (n(q)n eS (z) 


n(-q) ) = ( n(q)n eS (z)V. 


M~<l) 


z — S4 efI (z) / r \ n z — tt cS (z) 

i(q)n eS (z)Q n = (n(q)Ll eff (;z)n(-q)) r (n(q)n(-q)) 1 7 1 (n(q) ^ _ ^ cS ^ n{-q)j 


n(q)r2 ett (z)Q- 


' z - Q n fl eS (z)Q n 


Q n n ett (z)n{- q) ) (n(q)n(-q)) r ( n( q) 


z — n eff ( 2 ) 


*(-q) 


(58) 


r 


The important part of the first term on the right- To arrive at the penultimate line of Eq. (1471) we utilized 
hand-side of the last equality sign in Eq. (1551) is Eq. (1551) and to arrive at the last line we used the fact 


the matrix element of the effective evolution operator, that the local steady-state average followed by averaging 
(n(q)f2 eff (,j)n(—q)) , which can be expressed in terms of over particles’ positions is equivalent to the full steady- 
the frequency matrix 7i(q; z ), state average, 

(n(q)fr ff (*)n(-q)) r = ~q 2 NH(q ; z). (59) 

The frequency matrix is given by the following expression 

n(q-,z)= (60) 


«-U r = <•••>• 


(62) 


q ■ (^2i,j (m, sa < f *)lss (fj)lss) e 

Wo (z + Tp 1 ) 


-iq-(ri-rj) 


q 


In Eq. (l60l) q is a unit vector, q = q/q. In turn, the non¬ 
trivial part of the frequency matrix is represented by the 
function W||(g) introduced in Sec. [Vl which quantifies 
correlations of velocities of individual particles, 


We note that the penultimate line of Eq. (S3 shows 
that W||(g) quantifies fluctuations around the force bal¬ 
ance condition 


m 2 0 


w 

i 

W 


q- (E « f * f Ass - (^)iss (W J ) q = M(q; z)t 


The important part of the second term at the right- 
hand-side of the last equality s ign in Eq. (1551) can be 
expressed in terms of reducible |36t 1371 1 memory matrix 




q- (E + F >) « f 4>is S + q = 


q-(E( f i + F 0( f i+ F j)e" iq ' (ri “ rj) }q = W||( 5 ) \ Z ~ Q n n° s {z)Q ri 


n(q)n ett (z)Q 


Q n n etl (z)n(- q) ) = 


h3 


q 2 NM{q\z). (63) 

(61) The memory matrix is given by the following expression 


M(q; z) = (iV£ 0 4 (z + r- 1 ) 2 ) ' q . (e e- iq ' ri ~ (V lss <fj)J 

\ iJ 


h Vj + Ff ] Qr 


Z - Q 


Qeti/..\Q O' 71 E ~ ((fifro)lBH (f,) lss (f m )iJe-E q 
1 n l,m / 


(64) 


We can now rewrite the Laplace transform of the in¬ 
termediate scattering function in terms of the frequency 


and memory matrix, 


F(q;z) 


_ M _ 

z + q 2 (' H{q ; z) - M{q; z)) /S(q) 


(65) 
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where S(q) is the steady-state structure factor, 

S(q) = <n(q)n(—q)) r = (n(q)n(-q)) . (66) 

The second equality sign in Eq. (l66l) follows from the fact 
that for self-propulsion-independent quantities averaging 
over particles’ positions is equivalent to averaging over 
the full steady-state distribution of positions and self¬ 
propulsions. 

Eq. (1651) is the active matter equivalent of the mem¬ 
ory function representation of the intermediate scatter¬ 
ing function of a Brownian colloidal suspension derived 
by Ackerson [35?]. As pointed out by Cichocki and Hess 


[36| and later elaborated by Kawasaki [?37j], the memory 
matrix that enters into Ackerson’s formula is not fully 
irreducible. Both Cichocki and Hess, and Kawasaki ar¬ 
gued that approximations should be applied to the fully 
irreducible memory matrix. Therefore, we will follow the 
analysis presented in Refs. [36, 37] and derive an irre¬ 
ducible memory matrix. 

We define the irreducible evolution operator fi lrr (z), 
n iII (z) = Q n n eS (z)Q n ^sn iII (z) (67) 

where the subtraction term <5fi lrr (z) |?38j reads 


cKl irr (z) 


Qn (a^o 4 (z + T- 1 ) 2 ) ' £ V, ■ ((f,f,) lss - <f)) lss <f,)J e* qrj \ • q (n(q; z))~ 1 

ij / r 

q • (E e_iqrfe (M*. ^ ftU • (V« - PFf) Qn 

\ k,l 


( 68 ) 


Next, we define the irreducible memory matrix A4 lrl (q; z), which is given by the expression analogous to Eq. (1641) but 
with the projected evolution operator Q n f2 eff (z)<2„ replaced by irreducible evolution operator H lrr (z), 

M^(q; z) = (iV£ 0 4 (z + r- 1 ) 2 ) _1 q • (E e_iq " ((W\ ss - (fi)iss (fj)J ■ 

\ iJ 

x [—Vj + Ff ] Q n —^—Q n J2 V, • «f)f m ) lss - (f*) lss (U lss ) e iq ' Fm \ ■ q. (69) 

^ ' l,m / 


Finally, we use an identity similar to Eq. 


Z- Q n n eS (z)Qn z-n^{z) v ' 

I_^_ <TQeff / \_^_ 

+ z-a^(z) ou [z> z-Q n neZ(z)Q n ' 

and we derive the following relation between A i(q; z) and 
M ilr (q;z), 


M(q; z) = A4 irr (g; z) - M m {q ; z)H~ l {q] z)M(q; z). 

(71) 

Combining Eqs. (|65l) and (17T1) we arrive at the follow¬ 
ing representation of the intermediate scattering function 
in terms of the irreducible memory matrix, 


F{q-z) 


S{q) 

q 2 H(q;z)/S(q) 

1 +M" r (q-,z)/'H(q-,z) 


(72) 


Eq. (l72l) is the active matter equivalent of the mem¬ 
ory function representation of the intermediate scatter¬ 
ing function of a Brownian colloidal suspension derived 
by Cichocki and Hess [36j. The latter equation was the 


starting point of the derivation of the mode-coupl ing the¬ 
ory for the glassy dynamics of colloidal systems [39[ . 

In spite of the formal similarity between Eq. m and 
the irreducible memory function representation of the 
intermediate scattering function of a Brownian system, 
these equations are only equivalent if we take the limit 
z —> 0 in the frequency matrix term. In fact, as we have 
already stated a couple of times, active system performs 
interacting persistent Brownian motion and, therefore, 
its short time dynamics is ballistic rather than diffusive 
(we emphasize that, in this context, the term ballistic 
does not imply inertia, which is irrelevant for our over¬ 
damped system). The ballistic character of the short- 
time dynamics can be clearly seen if Eq. m is re-written 
in the following way. First, we rewrite Eq. (1721) by intro¬ 
ducing the irreducible memory function M lrl (q ; z), 

(z + r" 1 + M lrr (g; z)) (. zF(q ; z) - F(q; t = 0)) 

= - (u\\(q)q 2 /S(q)) F(q-,z). (73) 

where the memory function reads 

M hr (g;z) = (z + r- 1 ) 2 M irr (g;z)/o; n (q). (74) 
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Next, we rewrite Eq. (1731) in the time domain, 

d?F(q; t) + T~ 1 d t F(q- 1 1) + f (<K *) = 

- f dt'M ilr (q-t-t')d tl F(q-t'), (75) 

Jo 

where M 1IT (q;t) is the inverse Laplace transform of 
M lrr (g; z). We note that the resulting equation, Eq. (1751) . 
has the same form as the memory function equation for 
the intermediate scattering function of a Newtonian sys¬ 
tem. 

Let us discuss the meaning of the different terms in 
Eq. (1751) . starting with the three terms at its left-hand- 
side. The presence of the second time derivative implies 
ballistic dynamics at short times. For our system of self- 
propelled particles the origin of the ballistic behavior is 
the persistence of the microscopic motion. The first time 
derivative term describes the relaxation of the ballistic 
motion due to the evolution of the self-propulsion. Since 
we neglected the influence of the interactions on the evo¬ 
lution of the self-propulsion (recall Eq. (13910 . the relax¬ 
ation time of the self-propulsion is unchanged and equal 
to T p . The third term at the left-hand-side of Eq. m 
describes the collective random motion of the particles 
on time scale longer than the persistence time of the self¬ 
propulsion. On such a time scale, the second derivative 
term in Eq. can be neglected. Then (recall that we 
are neglecting the right-hand side term for a moment), 
the intermediate scattering function F(q;t) evolves dif¬ 
fusively, 

F{q-1) oc exp ^ • (76) 

and W|| (q)T p /S(q) plays the role of the short-time collec¬ 
tive diffusion coefficient. This simple exponential time 
dependence of density correlations is modified by the 


presence of the memory function term at the right-hand- 
side of Eq. m- The memory function term describes 
an internal time-delayed friction generated by the inter¬ 
particle interactions. This term leads to slow and glassy 
dynamics through a feedback mechanism implicit in the 
mode-coupling approximation (see discussion after Eq. 
»• 

The non-equilibrium nature of the active system man¬ 
ifests itself in Eq. (1751) through the presence of W||(g), 
which quantifies spatial correlations between velocities of 
different particles. In particular, the second time deriva¬ 
tive of the intermediate scattering function at t = 0 is 
expressed in terms of W||(g), 

dtF{q\t )| t=Q = F(q;t = 0) = -wy (q)q 2 , 

q _ (77) 

Eq. (1771) agrees with the short-time expansion discussed 
in Sec. 0 In the next section we will show that wy ( q ) 
also enters an approximate expression for the so-called 
vertex function. 

VII. MODE-COUPLING-LIKE 
APPROXIMATION 

In order to make Eq. (1751) useful we need an explicit 
expression for the memory function. Since the memory 
function contains all the non-trivial dynamic information 
about the system, there is little hope to derive an exact 
expression for it. Here we use the factorization approxi¬ 
mation, which is also at the heart of the mode-coupling 
theory of glassy dynamics, [l6j and thus we call our ap¬ 
proach a mode-coupling-like approximation. 

Specifically, to derive an approximate expression for 
the memory function we follow the steps of the derivation 
of the mode-coupling theory for systems evolving with 
Brownian dynamics [39(. The derivation consists of three 
steps. 

First, we project onto the subspace of density pairs, 


q . e- qr - «f t f,> lss - <f)) lss <f,)J • [-V, + Ff ] Q n 
\ ij 

« E q • ( E e_lqri - <5>1» (ff)lJ • [-V, + Ff] Q„n 2 (—qr, — q 2 ) 


x [(Qn7l 2 (qi,q 2 )Q„n 2 (-q 3 ,-q4)) r ] 1 (Qn«2(q3,q4) • 

I- 


(78) 


Here n 2 (qi,q 2 ) is the Fourier transform of the micro¬ 
scopic two-particle density, 

n2(qi,q2) = E e_lqiri_lq2 ' rm ’ ( 79 ) 

l,m 


and [(Q„n 2 (qi,q 2 )Qn« 2 (-q 3 ,-q 4 )) r ] 1 is the inverse 
of the correlation matrix of microscopic pair densities. 
We should emphasize the importance of using in Eq. 
(1751) only the parts of the microscopic pair density that 
are orthogonal to the microscopic density, i.e. using 
Qn^ 2 (qi,q 2 ) rather than n 2 (qi, q 2 )■ The presence of the 
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operator Q n is necessary for the existence of the inverse of 
the correlation matrix (Q„n 2 (qi, q 2 )Qn« 2 (-q 3 , -q 4 )) r 

M- 

Second, we factorize averages resulting from substitut¬ 
ing projection (1781) into the expression for the memory 


J 


function and at the same time replace the irreducible op¬ 
erator D lrr (z) by effective evolution operator fl eff (z). We 
should emphasize that this factorization has to be done 
in the time domain, 


CT~ l 

cr~ x 


Q ra n 2 (qi,q 2 ) (z - ^ lrr 0)) 1 Q„n 2 (—q 3 , -q 4 ) 


(80) 


n{ qi) (z - H 6 ®^)) 1 n(-q 3 ) 


CT 


i(q 2 ) (z - D eff (,z)) 1 n(-q 4 )) +{3o-4}. 


Here Cl 1 denotes the inverse Laplace transform and 
{3 4} means the preceding expression with labels 3 

and 4 interchanged. Consistently with Eq. (1801) we also 
factorize the steady-state correlation matrix of micro¬ 
scopic pair densities and for its inverse we get 


[(Qnn 2 (qi,q2)Q„n 2 (—q 3 , -q 4 )), 


i-i 


(81) 


(n(qi)n(-q 3 )) r 1 (n(q 2 )n(-q 4 )) r 1 + {3 -B- 4} . 

Third, we approximate the vertex functions. Due to 
the presence of the velocity correlations, this last step is 
somewhat more complex than the approximation used in 
the derivation of the standard mode-coupling theory [39j. 
We will explain it on the example of the left vertex, V;, 
which is given by the following formula 

Vz(q; qi,q 2 ) = 

£o" 2 q • (E e_iq ri ~ < f *>iss <fi) J ■ 

\ i,j 

x [—Vj + F° ff ] S„n 2 (- qi ,-q 2 )) r 


= c 2 q-(E e “ iqri ( f * + F »)( f i + F i)' 

\ i,j 

x [-Vj+F® ff ] Qnn 2 {~qi, -q 2 )) 


(82) 


Due to the presence of the projection operator Q n , the 
left vertex consists of two terms, 


V;(q; qi,q 2 ) = 

- ?o" 2 q- /E' 


(83) 




(f» + Fj) (fj + Fj) ■ 




x [Vj-n 2 (— qi , —q 2 )]) 


+ 2 q • ( E 1 




(fj + Fj) (fj + Fj) ■ 


1,3 


x [Vj-n(-q)]) (n(q)n(-q)) 1 (n(q)ra 2 (-qi, -q 2 )). 

We will start with the second term, which is a bit easier 
to analyze. With the help of convolution approximation 


1111 generalized to the stationary state of our active sys¬ 
tem, 


(n(q)n 2 (-qi,-q 2 )) = EE^«^ 

\ * j,k / 

« NS(q)S( qi )S(q 2 )6^ i+q2 j (84) 

we can rewrite the second term in the following form, 

£o“ 2 q-(E e “ iq ' ri (f, + F 8 )(f J + F J ). 

\ i,j 

x [Vj-n(-q)]) (n(q)n(—q)) 1 (rc(q)n 2 (-qi, -q 2 )) 

« *C 2 q • ^E ^ + F ‘) ( f i + F i) • q 

xS(q 1 )S(q 2 )6 cl 

,qi+q2 

= iTVcc||(<7)q • (qi + q 2 ) S(qi)S(q 2 )5q ,qi+q 21 (85) 

where in the last step we used the fact 
that for an isotropic system the tensor 

( f i + F *) ( f j +Fj)e -M1 ' (ri-r ^ has only com¬ 
ponents parallel and perpendicular to q, and that its 
parallel component is proportional to W||(q), see Eq. 

(S3. 

The first term at the right-hand-side of Eq. (1551) can 
be rewritten as follows, 


-&> 2 q ■ (E 1 




(fi + Fj) (fj + Fj) • (86) 


i,3 


x [Vj-n 2 (— qi) — q 2 )]) 

= -?o _2 q ■ ( E e _<q ' ri ( f < + F «) & + F ^)' 


i,3 


*qie iqi - r ^e iq2 ' ri +{1^-2} 


Now, we use an approximation which is a generalization 
of the convolution approximation to correlation functions 
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involving active particles’ velocities, 

£o” 2 (E e " iqr ‘ ( f *+ F >) ft+ F ^) eiqirj E eiq2 ' r ')« 

\ i,j l / 

lV£ 0 2 u;(q) • ((^f))" 1 •a;(q 1 )5( g2 )5 q 

,qi+q2- (87) 

Here w(q) is the tensorial version of W||(q), 

w (q) = (E & + F ‘) (fi +Fi)e- iq - (r *-^^ ,(88) 

and ((<5f<5f)) — 1 is the inverse of the correlation matrix of 
single-particle velocities, 

(sm) = 1^E & + F *) + F 0^ ■ (89) 

Note that lim^oo w(q) = (<5f<Sf) /£q. Using convolution 
approximation (1571) and symmetry properties of tensors 
cj(q) and (6f6f) we can finally write the first term (1551) 
in the following form 

-iN q- [qiW||(g)w||(qi)5'(fc 2 )/w||(oo) 

+q 2 W|| (g)w|| (q 2 )S{q 1 )/u\\ (oo)] <5 q . qi 

+q2 > (90) 

where W||(oo) = lim^oo W||(g) = 

(3-^Co) _1 (j2i (ft + F *) 2 )- 

Combining Eqs. (1551) and (1501) we get the following 
approximate expression for the left vertex, 


V;(q;qi,q 2 ) « -iNS(qi)S(q 2 )u\\{q) 

wii(gi) 


qqi 


- 1 + {1 O 2} 


^ w || (°o)S(qi) 

= iNpS(q 1 )S(q 2 )oj\\ (q) [q • qiC(gi) 4- q • q 2 C(g 2 )](91) 


where a new function C{q) reads 

wii(g) 


pC(q) = 1 - 


W|| (oo)5(9)' 


(92) 


We note that if correlations between particles’ velocities 
are neglected ( i.e . if W||(g) is set equal to 1 ), C(q) be¬ 
comes an analogue of the direct correlation function for 
our steady-state active system. 

The right vertex can be analyzed in the same way. We 
note that in the standard mode-coupling theory there is 
a close correspondence between the vertices that enter 
the expression for the memory function for the collec¬ 
tive intermediate scattering function and those that en¬ 
ter the expression for the memory function for the self¬ 
intermediate scattering function. We should emphasize 
that this correspondence is also present in our approach; 
compare Eqs. m and (IB27I) . Since the derivation 
of the vertices that enter the expression for the mem¬ 
ory function for the self-intermediate scattering function 


is more straightforward, the presence of this correspon¬ 
dence serves as an additional check for the derivation 
presented above. 

Combining the three steps and taking the thermody¬ 
namic limit we arrive at the following expression for the 
irreducible memory function, 

«■"(,;*) = (93) 

X (q- [qiC(gi) + q; 2 C(q 2 )}) 2 F( qi ;t)F(q 2 ;t). 

Expression (1551) resembles the corresponding expres¬ 
sion derived within the standard mode-coupling theory. 
It describes a dynamic feedback mechanism: the time- 
delayed internal friction arising due to interparticle inter¬ 
actions decays due to the relaxation of the two-particle 
density, which is included at the level of factorization 
approximation (1551) . Thus, slow decay of density fluctua¬ 
tions feeds back into slow decay of the irreducible mem¬ 
ory function. 

Eq. (1931) incorporates the non-equilibrium nature of 
the active system through the presence of the velocity 
correlations described through function wy ( 9 ). This func¬ 
tion sets the overall scale of the memory function. More 
importantly, it enters into the expression for new func¬ 
tion C(q), Eq. (1551) . and thus contributes to the vertices, 
which weight relative contributions of density fluctua¬ 
tions to the memory function. As we will see in the next 
section, it is the presence of the velocity correlations in 
the vertices that influences the location of the ergodicity 
breaking transition. 


VIII. ERGODICITY BREAKING TRANSITION 
AND LONG-TIME DYNAMICS CLOSE TO THE 
TRANSITION 

In general, one needs to solve the combined set of Eqs. 
m and (1931) numerically. However, as we show in this 
section, even without a numerical solution we can draw 
some general conclusions from these equations. 

First, let us assume that the system, as described by 
Eqs. (1751) and (1551) . undergoes an ergodicity breaking 
transition. At such a transition, the intermediate scat¬ 
tering function does not decay to zero. The long-time 
limit of this function, lim^oo F(q; t) = F(q) = f(q)S(q ), 
is the order parameter of the non-ergodic state. Follow¬ 
ing the procedure used in the standard mode-coupling 
theory [Hj we can derive a self-consistent equation for 
the normalized order parameter, f(q), 

= m(,) (94) 
where 771 ( 9 ) is given by the following equation 

™ 0 ?) = 7^1 J ^^ 2 <Hq - qi ~ c l2)S(q)S(qi)S(q 2 ) 
x (q • [qiC( 9 i) + q 2 C(g 2 )]) 2 f(qi)f(q 2 ). (95) 
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The self-consistent equation for the order parameter 
is very similar to the equation derived in the standard 
mode-coupling theory. The only difference is that the 
role of the function describing the effective interaction, 
which in the standard mode-coupling theory is played 
by the direct correlation function, is now played by the 
new function C(q) which involves both the steady-state 
structure factor, S(q), and the function describing spa¬ 
tial correlations of the particles’ velocities, Wj|(g). The 
presence of wy ( q ) in the vertex function means that the 
location of the ergodicity breaking transition is not ex¬ 
clusively determined by the local steady-state structure. 

Finally, to look at the long-time dynamics in the er- 
godic phase but close to the ergodicty breaking transi¬ 
tion, it is convenient to introduce a normalized correla¬ 
tor, (j>(q;t) = F(q-t)/S(q). The equation of motion, Eq. 
(iia, in the Laplace space and re-written in terms of the 
correlator (j){q ; z), has the following form, 

<t>(q;z) _ z + t - 1 + (a>|| (q)q 2 /S(q))m(q; z) 

1 -z<t>{q\z) U\\{q)q 2 /S(q) 

where m{q\ z) is the Laplace transform of the reduced 
memory function m(g; t), 

m (q-,t) = ^2 J d ^3 2 (5 ( q ~ qi ~ <i2) s {q) s (qi)S(q2) 

x (q- [qiC(gi)+ q 2 C(g 2 )]) 2 ^(q , i;t)</>(g2;^97) 


Near the ergodicity breaking transition, for small z the 
memory function m(q ; z ) becomes very large. Thus, for 
small z ( i.e . for long times), we can approximate Eq. 
(1M)| by the following one, 


<Kg; z) 

1 - Z(j){q\ z) 


m(q; z). 


(98) 


Eq. (1551) has the structure identical to that of the equa¬ 
tion of motion for the long-time dynamics near the er¬ 
godicity breaking transition described by the standard 
mode-coupling theory. This means that all analytical re¬ 
sults based of the standard mode-coupling theory carry 
over to the present theory for the dynamics of the active 
system. In particular, the only quantity that one needs to 
do in order to predict the so-called mode-coupl ing expo¬ 
nents is to calculate the exponent parameter A [161 ] . This 
parameter can be calculated from the solution of the self- 
consistent equation for the order parameter, Eqs. (19411951) 
at the ergodicity breaking transition. 

The reduced memory function m(q; z) differs from the 
corresponding quantity of the standard mode-coupling 
theory by the presence of the new function C(q ), which 
involves S(q) and W||(g). Again, due to the fact that 
the correlation function of particles’ velocities enters into 
m(q;t), the static structure factor does not completely 
determine the system’s dynamics. 


IX. DISCUSSION 


We presented here the details of the derivation of a re¬ 
cently proposed theory for the dynamics of dense ather- 
mal active systems. The theory, in a natural way, iden¬ 
tifies a new function that influences both short- and 
long-time dynamics. This new function quantifies spatial 
correlations of the particles’ velocities in the stationary 
state. The presence of the new function implies that the 
dynamics is not determined by the local structure only. 

The influence of the additional steady-state function on 
the dynamics allows us to describe an un-expected result 
obtained in computer simulations [2l|. We found that 
for a range of single-particle effective temperatures upon 
increasing departure from equilibrium (by increasing the 
relaxation time of the self-propulsion), the local struc¬ 
ture becomes monotonically more pronounced whereas 
the long-time dynamics first speeds up and then slows 
down. We showed in Ref. [2l[ that our theory, com¬ 
bined with steady-state correlation functions obtained 
from computer simulations of a model single-component, 
moderately supercooled system, is able to describe qual¬ 
itatively correctly the non-monotonic dependence of the 
long-time dynamics on the persistence time of the self¬ 
propulsion. 

Our theory shares the most important approximation, 
the factorization approximation, with the mode-coupling 
theory of glassy dynamics. We showed that if there is 
an ergodicity breaking transition, our theory predicts 
that the dynamics upon approaching this transition is 
qualitatively similar to that predicted by the standard 
mode-coupling theory close to the corresponding mode¬ 
coupling transition. 

In future work we would like to extend the present 
approach to binary mixtures in order to test the theory 
quantitatively against computer simulations, along the 
lines of our earlier test of the standard mode-coupling 
theory (42J. Of particular interest would be to check 
whether the present theory can predict the location of the 
dynamic crossover in an active system more accurately 
than the standard mode-coupling theory can predict the 
location of the crossover in an equilibrium (thermal) sys¬ 
tem. 

In addition, we would like to investigate the very fun¬ 
damental approximation of the theory concerning the ab¬ 
sence of the steady-state currents. In this context, we 
would like to develop and test a theory for effective tem¬ 
peratures in the active system. 

Finally, in real experimental active systems both self¬ 
propulsion and thermal Brownian forces are very often 
present. For this reason, we would like to develop a more 
general theory that would be able to account for the pres¬ 
ence of the Brownian noise. 
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Appendix A: Alternative approximation for 

Si ss SlQi SS 

We start by noting that due to Eqs. (1571) and (1551) the 
quantity we need to calculate is a specific matrix element 


of the inverse operator (z — Qi ss SQi ss ) 1 

^ Vi • y dt 1 ...dfjv (fi - (fi) lss ) (z - QissSQiss ) -1 

i 

X E ft - &>J P N^ 1. f l. UV, fiv) ■ V, 

3 

= E ' (& - (* - Qiss^Qlss )' 1 (fi - (f,) lss )) iss 

ij 

xPx(r lt ...,r N ) ■ V m . (Al) 

Following the spirit of the first Enskog approximation 
used in the kinetic theory we introduce the following ap¬ 
proximation, 


E Vi • ((fi - <fi) lss ) (z - QissHQiss)- 1 - (fi) J ) iss P%(n,...,r N ) • V m 

* E V ’ - (< f ^>lss - ( f *)l SS ■ [<fe - <fi>J (* - Qlss^Qlss) (fi - (fi) lss )) J" 1 

■ ((fif m ) lss - (fi) lss (f m ) lss ) P S N S ( r u ...,r N ) • V m (A2) 

Then we note that 

<(fi- - <fi> J 21-^Qiss (fi - <fi),J> lss = <(fi - (fi) lss )^fi) lss (A3) 

= [PSiTu...,^)]- 1 J dfr-.cMfi - = r- 1 (<f,fi) lss - (fi) lss (fi) lss ) - 2 DfSjlSap. 

Eqs. (1A2I1A3I) lead to the following approximate effective evolution operator, 

E Vi ' ((fifii) lss ~ (fi) lss (fi?) lss ) ' [( z — T p 1 ) ((fi?fi)l ss — (fi')lss (fi)lss) + ZDfSjlSapl 1 

• «fiU lss - <fi) lss (fin) lss ) • [-Ff + V,] . (A4) 


we limit the comments to the minimum. 

We start be rewriting the Laplace transform of the 
self-intermediate scattering function in terms of the ap¬ 
proximate evolution operator (1431) . 

FT [F s (q; t)] (z) = F s (q; z) = 

(n s { q) (z - fiy 1 n s (- q)(> = 

(n s (q) (z - S cff ( 2 )) 1 n s (-. (Bl) 

To derive the memory function representation of 
F( q\ z) we use the projection operator approach [Ifl 13(1 
|39(. We dehne a projection operator on the microscopic 
tagged particle density 

V s = n s (- q)) r (n s (q)... . (B2) 


It can be showed that expression (1A4I) reproduces cor¬ 
rectly the third time derivative of the scattering functions 
at t = 0. However, at present it seems that additional 
approximations are necessary in order to use tiHD as the 
starting point for the derivation of the memory function 
representation for the scattering functions. 

Appendix B: Tagged particle (tracer) density 
fluctuations 

The derivation of the approximate equation of mo¬ 
tion for the self-intermediate scattering function, which 
describes the time dependence of the tagged particle 
(tracer) density fluctuations follows the steps of the 
derivation presented in Sections IVTll VIII In this Appendix 
we present all the formulas but, in the interest of brevity, 
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Note that (n s (q)n s (—q)) r = 1 and thus we do not need 
to include a normalization constant in definition (1B2I) . 
Next, we use the identity 


1 

Z — fl eff (z) 


1 


z - n eff (z)Q s 

1 

+ z-n^{z)Q ; 


(B3) 


n eS (z)v s 


z - f l cS (z) ’ 


where Q s is the projection on the space orthogonal to 
that spanned by the microscopic tagged particle density, 

Q s =l-V s , (B4) 

to rewrite the Laplace transform of the time derivative 
of F s (q;t) in the following way 


CT[d t F s (q-,t)](z ) = ( n s (q)fl ett (z) 


z- fl eS (z) 


(-q)) = (n s (q L )ft ett (z)P ! 


1 


z - fl eff (z) 

1 


s (-q) 


+ ^n s (q)n cff (z)Q s —^pyn a (-q)^ = (n s (q)f2 cff (z)n s (-q)) r (n s (q)— ^-^(-q)^ 

1 


+ ( n s (q)f2 eff (z)Q. 


Q s fl eff (z)ns(-q) ) ( n s (q) 


z-Q s n eS (z)Q i 


z - fl eff (z) 


z - L2 eff (z) 

n a (~ q) 


(B5) 


The tagged particle frequency matrix, 'H s (q ; z), is defined 
in terms of the matrix element of the effective evolution 
operator, 

(n s (q)fi cff (z)n s (-q)) r = -q 2 H s (q; z), (B6) 


where 


The tagged particle memory matrix is defined as follows, 


^n s (q)fi off (z)Q s z - Q s n cS (z)Q s 2* neff (*) n »( **))_. “ 

q 2 M s {q-,z). (B8) 


H s (g;z) 


q- (( f i f i)i ss - ( f i)i ss ( f i)i ss ) r -q 

£o (* + t p _1 ) 


^ll(oo) 

Z + Tp 

(B7) 


Explicitly, the memory matrix is given by the following 
expression 


J 


Ma(q\z) = (Co (z + T p 1 ) 2 ) q-(j2 e iq ' ri (< f l f J>l SS - ( f l)lss &> J • 

\ j 

X t-V.+Ff] Q s g _ Qs ^ {z) Qs QsY^Vi- (m* - (f ; ) lss (fi)iJ e iq ' ri ^ -q (B9) 


To define the irreducible tagged particle memory matrix where the subtraction term <5f2| ff (z) reads 
we start by introducing an irreducible evolution operator 
for the tagged particle motion, 

ft r (z) = Q s n eS (z)Q s - 6Qf(z) (BIO) 


J 


™ = Q S (eO^P -1 ) 2 ) 1 Z] ■ (( f * f l)lss - ( f i)lss ( f l)lss) e iq ' ri ) q Wsivz ))- 1 

i / r 

q ■ e_iqri - ( f l)lss Wi J • (Vz - PFf) Qs 


(BH) 





















17 


Next, we define the tagged particle irreducible memory matrix Al 1 / r (g; z ), which is given by the expression analogous 
to Eq. (IB9I) but with the projected evolution operator Q s fi cff (z)Q,, replaced by irreducible evolution operator Qf T (z), 

z) = (Co (z + r- 1 ) 2 ) ' q • «fif,> lss - <fi) lss ft) J • 

\ 3 

(7)2* E V ‘- « f ' f i)lss - (ft)lss (fi) lss ) e iq ' ri \ • q. 

I ! Y 




(B12) 


Finally, we use an identity similar to Eq. 

1 1 


z-Q s tt eS (z)Q s z-W?{z) 


(B13) 


1 


-snf(z)- 


z-n™( z y““ s y ~ J z - Q s n eS (z)Q s ’ 

and we derive a relation between M s (q', z) and yV!( rr (g; z ), 

M s (q;z ) = M i y(q\z) - M i y{q;z)'H~ 1 {q;z)M s (q-,z). 

(BM) 

The memory function representation for the self¬ 
intermediate scattering functions reads, 


F s {q\z) = 


1 


z + 


q 2 n s (q;z) 


(B15) 


1+M i "(q;z)/H s (q;z) 


where Ml ir (q-,t ) is the inverse Laplace transform of 
M‘ rr (q; 2 ). We note that for evolution on the time scale 
longer than the persistence time the second time deriva¬ 
tive term in Eq. (IB 171) can be neglected. The result¬ 
ing equation becomes equivalent to that describing the 
self-intermediate scattering function of a thermal Brow¬ 
nian system with a short time self-diffusion coefficient 
W||(oo)r p . The last observation suggests that, in order 
to quantify the effect of the time-delayed friction (repre¬ 
sented by the memory function term) on the evolution of 
the active system, one should normalize the time scale by 
the term describing the short-time dynamics, W||(oo)r p . 

Finally, we note that Eq. (IB 1 71) leads to the follow¬ 
ing expression for the second time derivative of the self¬ 
intermediate scattering function at t = 0 , 


Again, we introduce the tagged particle (self) memory 
function, 

Ms T (q;z) = (z + t- 1 ) 2 M m (q;z) s /u || (oo), (B16) 

we rewrite Eq. (IB 151) in the time domain, 
d?F s (q;t) +T~ 1 d t F a (q-,t) = (B17) 

-W||(oo )q 2 F s (q\t)- f dt'M£ r (q\t - t')d v F s {q\t'), 


d 2 F s {q]t) | t=0 = -W||(oo )q 2 F s (q\t = 0) = -W||(oo)g 2 , 

(B18) 

which agrees with the exact result (l5ll) . 

To derive an approximate expression for tagged parti¬ 
cle memory function, we again follow the three steps [39j, 
with some modifications due to the difference between the 
tracer (particle number 1 ) and all other particles. 

First, we project onto the subspace of joint tagged and 
fluid particle density, 


q • ( E e_iqri (Miss - (fi) lss <f,-> J • [-V, + Ff] Q s 
« E q • (E e_Jq ' ri «Miss - (fi)iss <f,> J ■ [-V, + Ff] Q s n s2 (— qi , —q 2 ) 

qiv»q4 \ j 

x [(Q s n s2 (qi,q 2 )Q s n s2 (-q 3 ,-q 4 )) r ] _1 (Qin s2 (q 3 , q 4 ) 


r 


(B19) 


Here n S 2 (qijq 2 ) is the Fourier transform of the micro¬ 
scopic joint tagged and fluid particle density, 

n s 2 (qi,q 2 ) = E e -iqi ' ri-iq2 ' rm , (B 20 ) 

ra> 2 


and [(Q 5 n s 2 (qi,q 2 )Q s n s 2 (-q 3 ,-q 4 )) r ] 1 is the inverse 
of the correlation matrix of microscopic joint densi¬ 
ties. It should be noted that in Eq. (IB 191) we are 
using parts of the microscopic joint density that are 
orthogonal to the microscopic tagged particle density, 
2in, s 2 (qi,q 2 ). The presence of the operator Q s is nec- 
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essary for the existence of the inverse of the correlation 
matrix (Q s ra s 2 (qi,q 2 )Q s n s 2 (-q 3 ,-q 4 )) r [13- 

Second, projecting onto joint tagged-fluid particle den¬ 
sities leads to an expression involving a four-particle cor¬ 


relation function which is factorized. This factorization 
approximation is the main approximation involved in the 
present derivation, 


CT 


-1 


Q s n s2 ( qi,q 2 ) (z - 1 Q«n s2 (-q 3 , -q4) 


(B21) 


CT 


— 1 


(qi) (z - n eff (z)) rc s (-q 3 ) 


CT 


-1 


(q 2 ) (z - Cl eS (z)) n(-q 4 ) 


r 


Consistently with Eq. (IB21I) we also factorize the steady- be rewritten as follows 
state correlation matrix of microscopic joint densities and 
for its inverse we get 


[(Q s n s2 (qi,q 2 )Q 3 ra s2 (-q 3 ,-q 4 )) r 
<*qi,q 3 ( rl (q2)n(-q 4 ))7 1 • 


(B22) 


The final, third step, is concerned with the vertices. 
The left vertex reads, 

V s ;(q; qi,q 2 ) = 

£o -2 q ■ (E e_iqri (( f i f i>i ss - ( f i)is S &> J ' 
x [-V i + Ff] Q s n s 2 (- qi ,-q 2 )) r 

= £ 0 “ 2 q ■ (E e_iq ' ri ( f i + F i) ( f i + F i) • 

x [-V,- + Ff] Q s n s2 (- q 1; -da)) (B23) 

Due to the presence of the projection operator Q n , the 
left vertex consists of two terms, 


-£o _2 q ■ ( E e_iqri ( f i + F i) & + f j) ■ 
x [V-,n s2 (— qi,-q 2 )]) 

= -4” 2 q ■ (E e " iqri + F i) & + f j) • 

iqie iqi - ri ^i^e iq2 ' ri +iq 2 e iq2 ^e iqi ' ri (l-^i) 

l> l 

= -4“ 2 q ■ (E e_iqri ( f i + F i) to + F ^) • 

\ 3 

/qi(: lqi ' r; f) ( i E e lq2 ' ri - *(qi + q 2 )^ie iq2 ' ri 

Z 

+iq 2 e lq2 ' r V qi - ri ]> 

-* [w||(oo)5(g 2 )q ■ qi - W||(oo)q- (q 4 +q 2 ) 

+ w || (<? 2 )q • q 2 ] ^q,qi+q2i 


(B25) 


where in the last step the following approximation was 
(B24) used 


e _ * q ri (fi + Fr) (fi + Fi) e iqi ' ri ]T < 


<(f 1 + F 1 )(f 1 +F 1 ))(e 


-i(q-qi)-ri p iq 2 -rj 


V s i(q; qi,q 2 ) = 

- Co _2 q• (E e_iq ' ri to + F i)to + F i) 
x [V j n s2 (-q 1 ,-q 2 )]> 

+ C 0 ” 2 q• (E e " ,q ' ri + F i) ( f i + F i) 

\ i 

x [Vj-n s (-q)]) (n s (q)n s2 (-qi,-q 2 )}. 

The second term can be expressed in terms of u\\ (oo) and approximate expression for the left vertex, 

V s z(q;qi,q 2 ) « -iS(q 2 )u>\\{oo) 
wii (q 2 


>*q2-r I 


>(B26) 


= W||(oo)5(g 2 ) 

Combining Eqs. (IB25I) and (IB25I) we get the following 


the steady-state structure factor, 


4- 2 q-/^e- lq ' ri (f 1 + F 1 )(f J +F J )- 

\ 3 

x [Vj 7 i s (—q)]) (n(q)n s 2 (—qi, —q 2 )) 

= iw||(cx))q- (qi +q 2 ) (S(q 2 ) - l)5 q , qi+q2 . 


x q • q 2 


- 1 


v W||(oo )S(q 2 ) 

= ipS(q 2 )u!\\ (q)q • q 2 C(g 2 ) 


(B27) 


where function C{q) = (1 — 1 /S(q))/p, Eq. (|92l) of the 
main text. The right vertex can be analyzed in the same 
The first term at the right-hand-side of Eq. (IB24I) can way. 
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Combining the three steps and taking the thermody¬ 
namic limit we arrive at the following expression for the 
irreducible memory function for the tagged particle mo¬ 
tion 

M s r i>T,t) = PW||(oo) J d(q — q, — q 2 ) 

x (q- q. 2 C(q 2 )) 2 F s ( qi ;t)F(q 2 ;t). (B28) 


As we mentioned in the main text of the article, the 
relation between the approximate expressions for the ir¬ 
reducible memory functions for the collective and tagged 
particle motion, Eqs. (1M1) and (IB28I) . is the same as in 
the standard mode-coupling theory. 
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